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NASA TECHNICAL PAPER 


RECOVERY OF PYROSHOCK DATA FROM DISTORTED 
ACCELERATION RECORDS 

I. INTRODUCTION 


High frequency shocks, such as those generated by pyrotechnic devices, often produce acceleration 
records (time histories) that are distorted. If the degree of distortion is great, the data is usually discarded 
as “bad” data. However, if the distortion is slight, the data is usually “cleaned up” by questionable, make- 
shift means. Often the “clean up” procedure introduces as much error as previously existed in the data. 
The purpose of this report is to outline techniques for “cleaning up” data so that a true signal is obtained 
and so that the “clean up” techniques will be reproducible in any laboratory. 

The majority of ordnance shock data is distorted by either or both of the following: (1) Base line 
shift, or (2) accelerometer resonance. In the following sections, characteristics, correction methodology, 
and examples of data will be presented. 

II. BASE LINE SHIFT 
A. Characteristics 

Base line shift occurs with piezoelectric accelerometers when shock levels approach 100,000 g’s 
(this may occur at much lower levels depending upon the make and model of accelerometer). This effect 
usually appears as a signal envelope axis of symmetry shift, that is, the time history envelope is symmetric 
about an axis other than the zero base line. Also, this effect may appear as a signal superimposed upon 
an exponentially decaying pulse or upon a sine wave. Combinations of these two shift varieties also occur 
frequently. 

It is believed by experts that saturation of charge amplifiers and/or temporary accelerometer crystal 
repolarization are responsible for inducing base line shifts. 


B. Correction Methodology 


Donald Baker Moore of Explosive Technology, Inc. (P.O. Box KK, Fairfield, California 94533), 
has begun experimentation in identifying and correcting distorted data. Accelerometer output (tune 
histories) (Fig. 1) may be integrated once to obtain velocity time histories (Fig. 1) and twice to obtain 
displacement time histories (Fig. 3). If a residual velocity of great magnitude or a continuously increasing 
displacement exists after the shock, a base line shift is assumed. The shift may not always be visible in the 
acceleration time history. Figures 1, 2, and 3 represent a distorted or total signal and its derivatives. Figures 
4, 5, and 6 represent a true signal and its derivatives. Notice that the velocity and displacement both go to 
zero after the shock. 



flCCE 


I NJ I 


-4000 

-B000 

-8000 

10000 


i 


1 2 3 

TIME -mi ] 1 i seconds 


TXF^l^ XIME — 1.2! 


-rxMEr — 


Figure 1 . Acceleration time history. 
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The distortion source can be assumed to be an erroneous bias pulse of the following form; 


ABias = ^oexPH/To) 


If the bias is a pure direct current shift, Tq ->■ and ag^^g = = a symmetric shift. Suitable values 

may be selected for the parameters, thus forcing the velocity and displacement values toward reasonable 
levels. Often the bias terms are so obvious that the constants Aq and Tq may be chosen by inspection. 

When this is not possible, mathematical averaging techniques may be employed to derive the constants. 
The total signal may be written as follows: 

^TOTAL^^TRUE + ^BIAS a jguE = ^TOTAL " ^BIAS . 


Any true pyroshock pulse may be expressed as follows: 


^TRUE “ 2. [exp(-aj,j t) sin(bjvj t)] . 

N=1 


Once the bias term, agj^g, has been derived, the true signal may be easily extracted. The total signal 
with bias may be written as follows: 


^TOTAL 


OO 


E 

N=1 


texpC-a^ t) sin(bjsj t)] 


+ Aq exp(-t/TQ) . 


C. Examples 

There are five commonly occurring cases of base line shift (including the case when agj^g = 0). 
Case No. 1 : A true signal, no adjustments are necessary (Fig. 7). 


^BIAS ~^o~^ 

^TRUE^^TOTAL • 

Case No. 2: A simple direct current shift where Tq ->■ <» and Aq = N (Fig. 8). 
^BIAS " N 


5 



aTRUE"^TOTAL~N . 

Case No. 3: A combined or complex shift, multiple adjustments are necessary. At t = 0, agj^§ = N, or 
Aq = N. By inspection the bias contains two exponentials, one forcing the signal toward zero from 
above, and one forcing the signal to zero from below. A computer curve fill will determine TqJ and 

To 2 (Fig- 9). 

^TRUE “ ^TOTAL ~ l®xp(-t/lQ j ) -t- exp(-t/Tp2)l • 

Case No. 4: By inspection it may be seen that the signal is stacked upon a sine wave. There is no direct 
current shift (Fig. 10). 

^TRUE ^ ^TOTAL “ sin(-t/TQ> . 

Case No. 5 : By inspection it may be seen that the signal is direct current shifted and stacked upon a sine 
wave. Combine case No. 4 and case No. 2 (Fig. 11). 


^TRUE “ ^TOTAL “ " sin(-t/TQ) . 
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Figure 10. Sine distorted acceleration time history. 


3 r 



-1 

0 10 20 30 40 

TIME 


Figure 1 1. Direct current-sine distorted acceleration time history. 
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III. ACCELEROMETER RESONANCE 


A. Characteristics 


Any environment with significantly high g levels at high frequency will excite an accelerometer into 
producing a resonant response. This response may appear as a high frequency wiggle superimposed upon 
the lower frequency real signal, or it may be completely hidden in the real signal. In reality, all frequency 
components of the acceleration signal above one-tenth of the resonant frequency of the accelerometer will 
be amplified by the accelerometer response function. The accelerometer can be mathematically approxi- 
mated by a very lightly damped simple harmonic oscillator. The accelerometer response function 
(accelerometer frequency transfer function) is prescribed if the damping factor and resonant frequency 
are known. The transfer function may be written as follows: 


"^ACCEL(“) “ 


1 


1 + 


2 i { CO 
COjs, 




1^ACCEL(")| " ^/^ACCEl^l ■ ''ACCELt‘^1 



The total signal is, in fact, a convolution of the real signal (true signal) with the accelerometer 
impulse response function (transfer function) in the time domain. 


B. Correction Methodology 

The major problem in correcting data is determining the accelerometer transfer function. The 
damping factor and resonant frequency of any given accelerometer may be determined by testing and 
analyzing the resulting data. Also, many accelerometer manufacturers provide characteristics for their 
instruments. Once the accelerometer frequency response transfer function is derived, the remaining process 
is straight forward. 

First, deconvolve the total signal by Fourier transforming the time signal into the frequency domain 
as follows: 


'^TOTAL (w) - F {^toTAl} 

Next, divide the total signal by the accelerometer transfer function to yield the true signal as follows: 


^TRUE “ 


^TOTAL 

Yaccel 


Finally, use the inverse Fourier transform to convert the true signal back into the time domain as 

follows: 


^TRUE - F ^ {yjruE (")} 
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~ ^ACCEL “ I ^ACCEL^^) I “ ‘^A ® ‘^A* 

(See Fig. 13.) 

By examining the total signal (Fig. 14), one may observe that it may be described as follows: 

YtotaL^^^ — ^TOTAL ~ ^A e ^ sin (0,995 w^t — 14.5) 

-0.005 CO At 

+ 0.6204 e ^ sin(0.5 co^t - 7.228] 

NOTE: For convenience coj,j of the system was described as a function of co^. 

Deconvolve the signal by transforming into the frequency domain as follows: 

^TOTAL^^*^^ ^ ^ i^TOTAL^^^H 

Next, divide the total signal by the accelerometer signal as follows: 
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YaCCEU'^A^' > "A Y;^CCEL<“/"A> 







This should be the true signal. Finally, transform back into the time domain to Obtain the “true” 
time signal as follows: 


^TRUE^t) - axRUE ” I^TRUE^^^i 


= 0.5 co^ e 


-0.005 co^t 


sin 0.5 w^t 


(See Fig. 17.) 


IV. SUMMARY 


In the past much usable data has been discarded or erroneously cleaned up. Now, it is possible to 
systematically and correctly clean up data so that the clean up procedures are reproducible. 

The two major sources of data distortion are accelerometer resonance and base line shift. These 
distortions may or may not be readily observable. In theory, they may be identified and removed yielding 
true data. 


V. CONCLUSIONS 


In theory distorted data may be readily corrected upon investigation by removing bias terms or 
accelerometer transfer functions or both. These clean up procedures remain to be proven in the laboratory 
or applied sense. 
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